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SECTION  1 


INTRODUCTION 


The  recent  introduction  of  a  damage  tolerant  design  requirement  for 
critical  structural  components  In  U.S.  Air  Force  gas  turbine  engines  has 
created  a  need  for  accurate  crack  growth  predictions  [1].  The  Engine 
Structural  Integrity  Program  (ENSIP)  requires  that  a  pre-existing  flaw  be 
assumed  present  In  a  component  and  that  fracture  mechanics  analysis  and 
testing  demonstrate  that  such  flaw  will  not  grow  to  a  catastrophic  size 
during  the  lifetime  of  the  engine.  To  meet  this  requirement,  crack  growth 
rates  must  be  predicted  under  typical  mission  spectra.  For  components  such 
as  turbine  disks  and  spacers  In  advanced  fighter  or  attack  aircraft,  the 
loading  spectrum  will  generally  consist  of  large  amplitude,  low  frequency 
fatigue  cycles  with  Interspersed  hold  times.  The  hold  times  have  durations 
of  tens  of  seconds  with  amplitudes  anywhere  from  the  maximum  of  the  fatigue 
cycles  to  less  than  half  that  amount,  depending  on  the  type  of  aircraft  and 
the  mission. 

Concurrent  with  a  design  philosophy  change  from  low  cycle  fatigue 
Initiation  to  a  fracture  mechanics  damage  tolerant  approach,  engine 
operating  temperatures  have  Increased  steadily.  This  has  led  to  the 
development  of  higher  temperature  nickel-base  superalloys  and  a  tendency  to 
use  these  alloys  under  conditions  where  time-dependent  material  behavior, 
becomes  Important.  Thus,  the  effect  of  hold  times  In  a  design  spectrum 
which  results  In  crack  growth  during  the  hold  times,  has  to  be  considered 
as  part  of  the  life  prediction  methodology.  Although  there  has  been 
substantial  effort  devoted  to  the  study  of  elevated  temperature  sustained 
load  crack  growth  In  numerous  materials  [2-6],  there  has  been  little 


research  directed  at  crack  growth  under  mission  spectra,  particularly  when 
such  spectra  involve  both  sustained  loads  and  fatigue  cycles  with  sustained 
loads  at  different  levels  than  the  maximum  load  amplitudes  of  the  fatigue 
cycles. 

In  a  previous  investigation,  a  linear  cumulative  damage  model  was 
applied  for  the  prediction  of  crack  growth  rates  under  simple  spectrum 
loading  consisting  of  fatigue  cycles  and  hold  times  [7].  For  Inconel  718 
at  649‘*C,  a  material  which  exhibits  a  significant  amount  of  time-dependent 
behavior,  it  was  found  that  the  linear  cumulative  damage  model  worked  well 
when  hold  times  were  at  the  same  level  as  the  maximum  load  of  the  fatigue 
cycle.  Here,  the  algebraic  sumnatlon  of  the  fatigue  crack  growth  rate  as 
determined  from  constant  load  amplitude  fatigue  tests  and  the  sustained 
load  crack  growth  as  determined  from  constant  load  crack  growth  tests 
provided  excellent  correlation  with  experiments  where  the  two  types  of 
loading  were  combined.  However,  when  the  sustained  load  level  was  dropped 
to  75  percent  of  the  fatigue  maximum,  the  linear  summation  model  broke  down 
completely.  The  net  growth  observed  experimentally  appeared  to  be  due  only 
to  the  fatigue  cycles,  and  the  linear  summation  model  grossly  overpredicted 
the  growth  rate  under  the  simple  spectrum. 

This  investigation  was  undertaken  to  evaluate  the  range  of  validity  of 
a  linear  summation  model  and  to  determine  the  contribution  of  hold-time 
loading  in  a  simple  spectrum  consisting  of  combinations  of  fatigue  cycles 
with  hold  times  of  various  duration  and  amplitude.  A  simple  analytical 
model  is  developed  based  on  the  concept  that  the  sustained  load  crack 
growth  is  retarded  by  the  application  of  the  overload  fatigue  cycles. 


SECTION  2 


EXPERIMENTS 

A  series  of  crack-growth-rate  experiments  was  conducted  using  the 
simple  load  spectrum  depicted  in  Figure  1.  The  spectrum  consists  of  a  single 
fatigue  cycle  of  frequency  IHz,  load  ratio  R=0.1  and  maximum  stress 
intensity  factor  of  either  40  or  50  MPa*m^.  Hold  times,  t^^,  of 
durations  varying  from  a  few  seconds  to  several  hundred  seconds  were 
applied  at  amplitudes  given  by  K  =  where  R^^  is  the  ratio  of  the 

amplitude  of  the  hold-time  to  the  maximum  amplitude  of  the  fatigue  cycle. 
Values  of  R^^  of  1.0,  0.9,  0.8  and  0.5  were  used  which  correspond  to 
overload  ratios  of  1.0,  1.11,  1.25  and  2.0,  respectively,  where  the 
overload  ratio  is  the  ratio  or  in  the  "overload"  fatigue  cycle  to  the 
baseline  K  of  the  hold-time  load.  All  tests  were  conducted  on  CT  specimens 
of  width  W  *  40  mm  and  thickness  b  *  10  mm  under  computer  controlled 
constant  conditions.  The  test  procedure  is  described  in  detc’il  in 
Reference  8.  Crack  lengths  were  determined  from  compliance  readings  taken 
with  the  aid  of  an  extensometer  during  the  unloading  portion  of  the  fatigue 
cycle. 

The  material  used  in  this  investigation  was  Inconel  718  at  649°C 
having  a  standard  heat  treatment  given  in  Table  1.  The  material  is 
identical  to  that  used  in  several  previous  investigations  where  it  has  been 
shown  that  it  exhibits  significant  amounts  of  time  dependent  behavior 
[7,8].  Several  test  conditions  were  evaluated  on  a  single  test  specimen 
under  constant  K  control.  For  each  condition,  the  crack  was  grown  until  it 
could  be  unequivocally  determined  that  the  crack  was  growing  at  a  constant 
rate.  This  was  accomplished  by  continuously  monitoring  average  growth  rate 
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using  a  least  squares  linear  regression  fit  to  the  crack  length  versus 
number  of  cycles  data.  When  no  statistically  significant  deviation  from  an 
average  growth  rate  was  obtained  from  additional  applied  cycles,  the  growth 
rate  was  defined  as  reaching  a  steady  state  value.  This  value  was  used  in 
the  plots  of  cycle  block  growth  rate  against  total  cycle  time  T  =  tj^+1  (see 
Figure  1).  In  most  cases,  the  total  crack  extension  used  to  establish  a 
steady  state  growth  rate  under  a  given  test  condition  was  less  than  1  mm 
and  was  typically  of  the  order  of  0.5  mm.  Thus,  over  30  test  conditions 
could  be  evaluated  on  a  single  specimen,  making  it  possible  to  duplicate 
and  verify  the  growth  rate  for  any  given  condition.  Growth  rates  in  a 
single  specimen  were  found  to  be  reproducible  to  an  accuracy  of  10  percent 
under  a  given  load  spectrum  using  K  as  the  crack  growth  rate  correlating 
parameter.  Both  the  reproducibility  of  the  data  and  the  achievement  of 
constant  growth  rates  under  constant  K  conditions  demonstrate  that  K  is  an 
adequate  correlating  parameter  for  this  material. 
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SECTION  3 


EXPERIMENTAL  RESULTS 

Crack  growth  rate  per  cycle  block,  da/dN,  was  plotted  as  a  function  of 
total  cycle  time,  T,  for  the  cases  of  =  40  and  =  50  MPa-m^  in  Figures  2 
and  3,  respectively.  Note  that  a  cycle  time  of  1  s.  corresponds  to  pure 
fatigue  cycling  at  1  Hz  with  no  hold  time.  The  symbols  in  each  figure 
represent  the  experimentally  determined  growth  rates  under  constant  K 
conditions  for  various  values  of  R^^.  The  three  solid  lines  in  each  figure 
represent,  from  top  to  bottom,  the  analytical  predictions  from  a  linear 
cumulative  damage  model  for  values  of  R^^  of  1.0,  0.9  and  0.8, 
respectively.  The  linear  model  numerically  adds  the  growth  rate  for  a 
single  cycle  of  1.0  Hz  and  the  contribution  from  sustained  load  crack 
growth,  a*t^»  where  a  is  the  growth  rate  per  unit  time  for  the  sustained 
load  crack  growth  at  a  stress  intensity  factor  value  K^.  Values  of  a  as  a 
function  of  K  were  determined  frcwn  a  parallel  investigation  using  the  same 
batch  of  material  [9].  Several  constant  load  tests  were  conducted  to 
obtain  crack  length  versus  time  data  which  were  reduced  to  a  versus  K  and 
fit  with  a  sigmoidal  equation  to  provide  numerical  values  of  a  for  any  K 
value. 

The  data  in  Figures  2  and  3  reconfirm  the  previous  observation  that  a 
linear  cumulative  damage  model  works  well  for  R^.^  =  1.0,  ie  when  the 
hold-time  load  is  at  the  maximum  value  of  the  fatigue  cycle.  For  values  of 
Rin  =  0.9,  the  linear  model  overpredicts  by  more  than  a  factor  of  2,  thus 
there  appears  to  be  some  retardation  of  the  sustained  load  crack  growth. 

For  values  of  R^^  =  0.8,  there  does  not  appear  to  be  any  contribution  of 
the  hold-time  loading  for  values  of  hold  times  up  to  200  s,  the  maximum 


used  in  this  investigation.  The  linear  model  grossly  overpredicts  the 
growth  rate  for  this  case.  The  same  observation  can  be  made  for  =  0.5. 
For  this  latter  case,  results  were  obtained  only  for  the  first  set  of 
experiments  where  =  40  MPa*m^.  There  was  no  reason  to  repeat  these  very 
time  consuming  experiments  for  =  50  MPa»m^. 


SECTION  4 


ANALYTICAL  MODEL 

A  simple  model  Is  developed  to  predict  the  crack  growth  rates  under 
the  simple  spectrum  of  Figure  1.  The  model  uses  a  linear  cumulative  damage 
concept  by  summing  the  contributions  of  the  single  fatigue  cycle  and  the 
hold-time  load.  The  fatigue  cycle  is  assumed  to  retard  the  sustained  load 
crack  growth  because  the  latter  occurs  within  an  enlarged  or  overload 
plastic  zone  due  to  the  fatigue  cycle.  This  retardation  concept,  similar 
to  that  used  by  Wheeler  [10],  is  illustrated  schematically  in  Figure  4. 

The  plastic  zones  due  to  the  sustained  load  of  stress  intensity  factor 
and  the  overload  fatigue  cycle  of  maximum  stress  intensity  have  radii 
denoted  by  rp  and  Fp  with  since  Kg/KjjjSl.  If  the  overload  is  applied 

when  the  crack  length  is  a^,  the  tips  of  the  plastic  zones  due  to  the 
baseline  (hold-time)  load  and  fatigue  overload  are  at  a^  and 
respectively.  For  a  crack  advance  of  Aa  from  a^^,  the  crack  will  be  growing 
in  an  overload  plastic  zone  until  Aa  =  a2  -  a^.  During  this  time,  the 
growth  rate  will  be  retarded.  The  amount  of  retardation  can  be  expressed 
by  an  effective  value  of  stress  intensity  factor,  Kgf^»  which  is  smaller 
than  the  applied  K.  Using  a  single-valued  function  of  a  versus  K  to 
represent  the  sustained  load  growth  rate,  a  value  of  *^eff  ^  ess  than  the 
applied  K  will  result  in  a  lower  growth  rate.  For  modeling  purpose,  "eff 
is  taken  in  the  form 

Keff  =  -  cxexp(-3Aa)j  (1) 


where  a  and  3  are  parameters  to  be  determined  and  Aa  is  the  amount  of  crack 
extension  after  application  of  the  overload  .  The  parameter  3  Is  chosen 


such  that  “  ^p”*^p*  steady  state  crack  growth  resumes 

when  the  crack  has  traversed  the  overload  plastic  zone.  This  is 
accomplished  by  setting 


3Aa  =  (2) 

when  Aa  =  resulting  in  a  value  of  which  Is  within  one  percent 
of  unity  as  shown  in  Figure  5.  The  plastic  zone  sizes  are  taken  to  be  the 
plane  stress  values  given  by 


:  =  1  fJk 
P  ^ 


(^) 


where  Is  the  uniaxial  tensile  yield  stress  of  the  material.  Combining 
(2)  and  (3)  provides  an  expression  for  g: 


khy"-i: 


where  the  overload  ratio  y  Is  defined  as 


Jii  =  ^ 


It  is  important  to  note  that  3  is  fully  defined  from  the  test  conditions 
and  material  properties  and  cannot  be  used  as  an  adjustable  parameter  to 
fit  experimental  data.  Thus,  eq.  (1)  has  only  one  parameter,  a,  which  can 


be  chosen  to  fit  data.  Referring  to  Figure  5,  it  can  be  seen  that  the 
value  of  a  determines  immediately  after  the  overload  application  and, 
therefore,  can  be  used  to  model  the  value  of  sustained  load  crack  growth 


SECTION  5 


DETERMINATION  OF  CONSTANTS 


Analytical  predictions  were  made  using  a  linear  cumulative  damage 

model  which  Incorporates  the  retarded  sustained  load  crack  growth  behavior 

given  by  eq.  (1).  Computations  of  crack  extension,  Aa,  were  made  by 

numerical  integration  of  da/dt  using  1  s  time  Increments  for  At.  The  crack 

growth  rate,  da/dt.  Is  expressed  In  the  form  of  a  modified  sigmoidal 

equation  fit  to  experimental  sustained  load  crack  growth  data  [9].  The 

stress  Intensity  In  the  sigmoidal  equation  Is  replaced  by  •^eff  from  eq. 

(1).  The  only  parameter  which  had  to  be  determined  was  a  and  this  was 

accomplished  by  trial  and  error  by  choosing  a  range  of  values  for  a. 

The  results  of  these  trials  are  shown  in  Figures  6  and  7  for  K  =  40  and 

m 

=  50  MPa*m*,  respectively,  and  for  the  case  of  ®  0.9.  The  plots  for 
a  *  0  represent  the  linear  summation  model  for  the  case  of  no  retardation. 

It  can  be  seen  that  values  of  a  =  0.15  and  a  *  0.27  fit  the  experimental 
data  for  the  two  cases  very  well,  respectively.  It  Is  noted  In  both  cases 
that  the  model  overpredicts  the  growth  rates  for  very  short  hold  times 
below  approximately  10  s. 

For  the  case  of  R^.^  =  0.8,  there  appeared  to  be  almost  complete 
retardation  for  hold  times  up  to  the  maximum  of  200  s  used  In  the 
experiments.  In  this  situation,  a  could  not  be  determined  accurately  as 

can  be  seen  in  Figure  8.  Complete  retardation  is  equivalent  to  a  zero  growth 

* 

rate  which  corresponds  to  a  value  of  '‘eff  •  ''  which  Is  the  threshold  value 
of  K  In  the  §  versus  K  representation  of  the  hold-time  data.  In  our 
modeling  of  the  data  using  a  modified  sigmoidal  equation,  K  =  24  MPa*m’. 
Values  of  a  =  0.25  and  a  -  0.4  represent  complete  retardation  for  the  cases 
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of  *  40  and  ®  50  MPa*tn*,  respectively.  To  represent  the  data  in 
Figure  8,  values  of  0.246  and  0.392,  respectively,  were  chosen. 

Calculations  were  made  of  the  overload  plastic  zone  size,  defined 
as  through  which  distance  a  crack  must  grow  under  sustained  load 

before  it  returns  to  its  steady  state  growth  rate.  For  the  four  cases 
corresponding  to  the  two  values  of  and  for  =  0.8  and  0.9  for  each 
case,  values  of  this  overload  zone  size  ranged  from  approximately  0.1  to 
0.5  mm.  Note  in  Figures  2  and  3  that  a11  of  the  experimentally  measured 
growth  rates  per  cycle  block  were  well  below  these  values.  Thus,  for  the 
simple  spectrum  used,  the  sustained  load  crack  growth  was  always  within  the 
overload  plastic  zone  and  occurred  at  a  reduced  rate  due  to  retardation. 

In  fact  most  of  the  sustained  load  crack  growth  occurs  for  small  values  of 
Aa/(rp-rp),  or  near  the  origin  of  the  x-axis  shown  in  Figure  5.  For  these 
experimental  conditions,  the  parameter  a  is  the  primary  factor  governing 
the  growth  rate  in  the  analytical  model.  The  parameter  3,  and  the  shape  of 
the  function  shown  in  Figure  5,  are  of  secondary  importance  for 

short  hold  times  under  approximately  100  s. 


SECTION  6 


DELAY  TIMES 

In  studies  of  fatigue  crack  growth  retardation  due  to  overloads,  the 

concept  of  a  number  of  delay  cycles  is  commonly  used.  Here,  we  introduce 

the  equivalent  concept  of  a  delay  time  for  sustained  load  crack  growth. 

Figure  9  presents  plots  of  the  calculated  value  of  Aa^,  the  growth  during  the 

sustained  portion  of  the  load,  as  a  function  of  the  hold  times,  t^,  for  the 

cases  where  =  0.9.  Also  shown  in  dashed  lines  is  the  crack  extension 
in 

that  would  occur  if  there  were  no  retardation.  The  delay  times,  or  the  net 
gain  in  time  because  of  the  retarded  growth  rates,  are  approximately  200  s. 
as  shown  in  Figure  9.  Delay  times  of  approximately  4000  and  8000  s  were 
calculated  for  the  cases  where  =  0.8  and  *  40  and  50  MPa*m*, 
respectively.  In  a  parallel  investigation  on  the  same  material,  delay 
times  of  approximately  4000  and  20000  s  were  experimentally  measured  for 
R^^  =  0.83  and  R^^  *  0.67,  respectively,  for  similar  values  of  [9]. 


SECTION  7 


DISCUSSION  AND  CONCLUSIONS 

For  the  simple  spectrum  consisting  of  a  fatigue  cycle  and  a  hold  time, 
for  the  material  used  In  this  Investigation,  the  hold  time  does  not 
contribute  to  crack  growth  except  when  its  amplitude  Is  close  to  that  of 
the  maximum  of  the  fatigue  cycle.  If  the  sustained  load  crack  growth  Is 
considered  to  be  retarded  by  the  single  cycle  overload,  it  appears  that 
overload  ratios  greater  than  approximately  1.2  severely  retard  growth  for 

3 

hold  times  as  long  as  10  s. 

A  single  parameter  retardation  model  for  sustained  load  crack  growth 
combined  with  a  linear  summation  model  which  adds  cyclic  and  (retarded) 
sustained  load  contributions  appears  to  be  adequate  for  predicting  crack 
growth  rates  under  the  simple  load  spectrum  used  in  this  investigation. 

The  model  correlated  data  covering  a  range  In  hold  times  up  to  200  s, 
although  the  model  was  slightly  conservative  for  short  hold  times  up  to 
approximately  10  s.  It  appears  that  in  addition  to  the  retardation  of  the 
sustained  load  crack  growth  due  to  a  fatigue  overload,  there  is  also  a 
slight  retardation  of  the  fatigue  crack  growth  due  to  the  sustained  load. 
This  latter  phenomenon  was  not  modeled  in  this  investigation. 

Most  of  the  retardation  modeled  In  these  experiments  covered  short 
holo  times  corresponding  to  small  crack  extensions  compared  to  the  size  of 
the  overload  plastic  zone.  Thus,  the  initial  growth  rates  directly  after 
an  overload  appear  to  be  modeled  well.  In  spite  of  the  fact  that  overloads 
are  repetitively  applied  before  the  crack  has  grown  even  a  small  distance 
through  the  overload  plastic  zone,  a  steady  state  growth  rate  was  achieved 
in  all  experiments  under  constant  K  conditions.  This  provides  further 
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evidence  for  the  use  of  K  as  a  correlating  parameter  for  sustained  load 
crack  growth  In  this  material. 
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TABLE  1 


HEAT  TREATMENT  FOR  ALLOY  718 


STEP  1:  Anneal  at  968®C  (1775‘’F)  for  1  Hour,  Then  Air  Cool  to  Temperature 

STEP  2:  Age  Harden  at  718“C  (1325®F)  for  8  Hours,  Then  Furnace  Cool  at 
56“C/Hr  (lOO^F/Hr)  to  621®C  (1150“F) 

STEP  3;  Age  Harden  at  621®C  (1150®C)  for  a  Total  Aging  Time  in  Step  2  and 
Step  3  of  18  Hours 

STEP  4:  Air  Cool  to  Room  Temperature 


Loading  spectrum  consisting  of  1  Hz  fatigue  cycle  and 
sustained  load. 
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L_MODEL  Rin=1.0;  Kmax=50  MPA^/ff 
L_MODEL  Rin=0.9;  Kmax=50  MPA^/M 
LMODEL  Rin=0.8;  Kmax=50  MPA^/R 
EXPERIM  Rin=1.0:  Kmax=50  MPAv^ 
EXPERIM  Rin=0.9;  Kmax=50  MPA>/M 
EXPERIM  Rin=0.8;  Kmax=50  MPA^/M 


Figure  3. 


TIME  (sec) 

Crack  growth  rate  data  plotted  against  total  cycle  time 
for  1C  =  50  MPa»ro®  as  a  function  of  different  R. 
values;*  Linear  model  data  is  also  given  in  the  figure. 


overload 


10'^  10'  10*-  10^ 

TIME  (sec) 

Figure  7.  Calculated  crack  growth  rate  for  different  values  of  a 
^or  l(  =  50  MPa»m*  and  R.  =  0.9.  Data  are  also  shown 
In  thFfigure. 


Do/Dn  ( 


Figure  8.  Calculated  crack  growth  rate  for  various  a  for  »  0.8. 
Data  are  also  shown  in  the  figure. 
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Rin=0.9;  ALPHA=0.27;  Kmax=50  MPAv^ 
Kmax=50  MPA^;  NO  RETARDATION 
Rin=0.9;  ALPHA*0.15:  Kmax-40  MPA^ 
Kmax.40  MPA^;  NO  RETARDATION 

/  / 


DELAY  TIME  (td) 


Figure  9.  Calculated  crack  advance  during  sustained  load  after  the 
application  of  the  overload  as  a  function  of  the  sustain 
load  time  for  R.  *  0.9  and  *  40  and  50  MPa*m*  as 
predicted  from  tne  retardatlorPmodel .  Predicted  delay 
times  (t.)  are  also  shown  in  the  figure. 


